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ABSTRACT 

Context. Lithium and beryllium are destroyed at different temperatures in stellar interiors. As such, their relative abundances offer 
excellent probes of the nature and extent of mixing processes within and below the convection zone. 

Aims. We determine Be abundances for a sample of eight solar twins for which Li abundances have previously been determined. The 
analyzed solar twins span a very wide range of age, 0.5 - 8.2 Gyr, which enables us to study secular evolution of Li and Be depletion. 
Methods. We gathered high-quality UVES/VL.T spectra and obtained Be abundances by spectral synthesis of the Be n 313 nm doublet. 
Results. The derived beryllium abundances exhibit no significant variation with age. The more fragile Li, however, exhibits a mono- 
tonically decreasing abundance with increasing age. Therefore, relatively shallow extra mixing below the convection zone is necessary 
to simultaneously account for the observed Li and Be behavior in the Sun and solar twins. 
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1. Introduction 


The light elements lithium and beryllium are fragile, meaning 
that they are destroyed at temperatures of about 2.5 ■ 10 6 K and 
3.5 ■ 10 6 K, respectively (through a and proton captures). The 
observed abundances of Li and Be thus provide constraints for 
the transport of material in stellar interiors. To decrease the Li 
and Be abundances, the material has to be transported to deeper 
and hotter regions within the star before returning to the surface. 
In the Sun, lithium destruction requires temperatures somewhat 
hotter than those achieved at the base of the convective zone ac¬ 
cording to standard models of stellar evolution. The fact that the 
observed photospheric Li abundance is some 150 times lower 
than the meteoritic value thus necessitates extra mixing below 
the convection zone. The mechanism (or mechanisms) for de¬ 
pleting the light elements is still debated. P ossib le processes 
leadin g to extra mixing incl ude rotation d Pinsonn eault et all 
19891) . internal gravity waves (ICharbonnel & Talonl 20051) . mi- 
croscopic diffusion and gravitational settling (Michaud et al. 
2004), and convective overshooting ( IXiong & DengH2007l) . Be¬ 
cause beryllium destruction requires greater temperatures than 
Li, its abundance serves to constrain the extent of this extra mix¬ 
ing. 


The solar photospheric Be abundance has been a source of 
contention. Early work suggested that Be was depleted in the 
solar photosphere compared to meteorites (Chmielewsk i et all 
i 19751 ), In solar-type stars, the Be abundances can only be esti¬ 
mated through the Ben doublet at 313 nm, a spectral region dif- 


* Based on observations obtained at the European Southern Observa¬ 
tory (ESO) Very Large Telescope (VLT) at Paranal Observatory, Chile 
(observing program 083.D-0871). 


ficult to analyze because of blends and uncertain atomic data. It 
has long been debated whether there is a substantial amount of 
missing UV opacity (e.g „ Magain 1987; Kurucz 199 2; Allende 
Prieto & Lambert 2000). lBalachandran & Belil d 19981) attempted 
to empirically calibrate the amount of missing continuous UV 
opacity by enforcing the same O abundances from the OH A- 
X lines near the Be doublet and the OH vibration-rotation lines 
in the infrared. They inferred a substantial amount of this miss¬ 
ing opacity and also that the solar Be abundance is the sam e as 
the meteoritic value within errors, a conclusion which lAsplundl 
( 2004 ) also reached using a more sophisticated 3D hydrodynam- 
ical model atmosphere. Without properly accounting for this ad¬ 
ditional continuous opacity, the Be abundance becomes under¬ 
estimated, which leads to erroneous conclusions whether there 
eve n is a solar Be depletion ; in fairness, we note, however, 
that lChmielewski et al.l d!975l) concluded that the uncertainty in 
their Be abundance (1.15 + 0.20 dex) was too large to advocate 
any substantial Be depletion. These suspicions about substantial 
missing UV continuous opacity were subsequently confirmed 
by Bell et al. (2001) using new calculations by the Iron Project 
for the bound-free opacity of Fe i. More recently, iTakeda et al.l 
(1201 ll) studied a sample of 118 solar analogs and suggested that 
Be depletion in the Sun could be significant. 


There have been several observational studies of Be in solar- 
type stars (Santos et aP 12004 [Boesgaard & Krugler Hollekl 
2009; lRandichll2010l: ITakeda et ail 20111) . but none focused on 
solar t wins, except for the qualitative work bv ITakeda & T aiitsu 
(2009) on three solar twins. The importance of solar twins is that 
as they have nearly solar mass and composition (iMelendez et al.l 
l2014bl) , their evolution is similar to that of the Sun. Hence, so¬ 
lar twins at different evolutionary stages in the main sequence 
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can be used as proxies of the Sun at different ages. In this work, 
we obtain Be abundances for solar twins in a broad age range 
to provide constraints on Be depletion during the main sequence 
and thus on the extent of any extra mixing below the convection 
zone. 


2. Observations and data reduction 

Spectra of eight solar twins and the Sun were obtained with 
the UVES spectrograph on the 8.2m UT2 Very Large Telescope 
at ESO Paranal, on 29-30 August 2009. We used the dichroic 
mode, obtaining simultaneous UV and optical coverage in two 
setups: i) with standard settings of 346 nm + 580 nm; ii) with the 
standard 346 nm and a nonestandard setting centered at 830 nm. 
We achieved a high S/N in the UV because the 346 nm setting 
(306 - 387 nm) was covered in both setups. The 580 nm stan¬ 
dard setting covered the optical (480 - 682 nm) region, and our 
830nm setting included the red region (642- 1020nm). The UV 
setup was used to obtain the Be abundance, the optical setup was 
used to obtain the stellar parameters and Li abundances. 

We achieved a resolving power R = 65,000 at 306 - 387 nm 
and R = 110,000 for 480 - 1020 nm. The solar spectrum was 
also obtained through the observation of the asteroid luno with 
the same spectrograph setup and served as the solar reference in 
our differential analysis. The spectral orders were extracted and 
wavelength calibrated using IRAF, with additional data process¬ 
ing performed with IDL. The reduced spectra have S /N ~ 100 
and 1000 near the Be and Li lines, respectively. 

There are no spectral regions free from lines in the near UV 
spectra of solar twins. Hence, we employed a dedicated contin¬ 
uum normalization technique, which is described in Ramirez et 
al. (2008). We took advantage of the superb continuum normal¬ 
ization of the solar spectrum reported by Kurucz et al. (1984), 
which is used as a reference. Each order of the UVES spectra 
was divided by its corresponding piece from the spectrum of Ku¬ 
rucz et al. after matching their spectral resolution, correcting for 
radial velocity offsets, and rebinning to a common wavelength 
sampling. In principle, the result should have been a smooth 
function corresponding to the shape of the continuum (the up¬ 
per “envelope”) of the UVES spectra. However, because of the 
finite S/N values and instrumental differences or small defects, 
this envelope had to be smoothed out using a 100-pixel wide 
(~0.16 nm) median filter. The UVES data were then divided by 
this envelope in each order. In essence, this procedure makes the 
UVES spectra inherit the continuum normalization of the spec¬ 
trum of Kurucz et al. 


3. Analysis of the Be abundances 


We used the doublet resonance lines of Be n at 313.0420 nm and 
313.1065 nm to determine Be abundances; these are the only 
Be lines available for observation from the ground in solar-type 
stars. As the lines are in the UV region and blended by different 
species of atoms and molecules, the abundance determination 

was carried out by means of spectral synthesis. _ 

Ou r initial line list was based on the list of AshweBetalJ 
(120051) and was checked with the list by Primasetal. (1 1 997l) ~ 
which is the list used bv lTakeda et alJ (1201 lit . In some cases we 
modified the log ^/-values to achieve a good agreement with the 
spectrum of the Sun and of a solar analog that is severely de¬ 
pleted in Be (Schirbel et al. 2015, in prep.). As the Be-poor solar 
analog does not show any Be lines, it was of great help to better 
constrain the lines blending the Be features in the Sun. We em¬ 
phasize that both the Sun and the Be-depleted solar analog were 




Fig. 1 . Comparison between the observed (points) and synthetic (lines) 
spectra of the solar twin HIP 102152 around the Be n lines. The best fit 
is the central line, the other lines show changes in the Be abundance by 
±0.05 and ±0.10 dex. 


observed with the same UVES setup as the solar twins and were 
reduced in the same way, allowing thus a reliable differential 
abundance analysis between the Sun and the solar twinsQln Fig. 
[Q we show the spectral synthesis for the solar twin HIP102152. 
The weaker Ben 313.1 nm line is less blended and offers more 
reliable abundances. Since the stronger 313.Onm line is heavily 
blended by CH and OH lines, abundances were determined for 
comparison, but not used in the final determination. 

For the spectral synthesis we used the synlli driver of the 
February 2014 version of the ID LTE code MOOG ( Sneden i 
119731) . which includes continuum scattering. We adopted A(Be) 
= 1 .38 dex as the standard solar Be abundance dAsplund et al.l 
2009). The model atmo spheres were interpolate d from the AT- 
LAS9 Kurucz’s grid (ICastelli & Kurucz] 120041) with the T e jf, 
l og#, [ Fe/Hl and micro t urbulen ce determined by iMonroe et al.l 
(12013h : IMelende z et al.l (l 2014al) . and Monroe et al. (201 5; in 
preparation). I Asplundl (12005b and PTakeda & Taii tsu (2009) con¬ 
cluded that the Be n lines are insensitive to non-LTE effects in 
the Sun. Considering the similarity in stellar parameters among 
the solar twins and the differential nature of the analysis, any dif¬ 
ferential NLTE corrections would probably be vanishingly small 
fe.g.. iMelendez et akll2012l) . hence no corrections were applied. 

To determine the macroturbulence line broadening, we first 
analyzed the line profiles of the Fe i 602.7050 nm, 609.3644 nm, 
615.1618 nm, 616.5360 nm, 670.5102 nm and Ni i 676.7772 nm 
lines in the Sun; the synthes es also included a rotational broaden - 
ing of v sin i — 1.9 km s~ 1 (iBninindll9841:ISaar & Ostenll 1 997h . 

1 Our modified log gf -values of the Ben lines are 0.6 dex lower 
than those recommended by iFuhr & Wiesel 1201 Oh . but we emphasize 
that our values are only valid for our internal differential analysis. In 
other words, they are only valid when using our set of reduced UVES 
spectra, the 20 14 version of MOOG, the set of solar abundances of 
lAsplund et all 120091) . and the adopted blends for the Be n region. More¬ 
over, our Be abundances are differential relative to the solar Be abun¬ 
dance of Asplund et al l l2009lf . 
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and the instrumental broadening. The macroturbulent velocity 
found for the Sun is V mac ro,0 = 3.6 km s 1 . For the solar twins, we 
estimate the macroturbulence following iMelendez et ali (12012|) 
(average of equations E.2 and E.3): V maC ro,star = V ma cro,o + (T e s - 
5777)/486. 

With the macroturbulence fixed, v sin i was estimated for the 
solar twins sample by fitting the profiles of the six lines men¬ 
tioned above, also including the instrumental broadening. Fi¬ 
nally, the best-fitting Be abundances were estimated using a^ 2 - 
procedure. The estimated macroturbulence, v sin i and Be abun¬ 
dances for the entire sample can be found in Table 1. 

We estimated the errors considering both observational and 
systematic uncertainties. The observational errors are due to un¬ 
certainties of the continuum placement and S/N (the synthetic 
spectra were shifted vertically within the allowed noise of the 
observed spectrum; the abundance variation due to this shift was 
adopted as the observational error). For the systematic errors 
we considered the errors in the stellar parameters. Both obser¬ 
vational and systematic errors were added in quadrature. 

Thanks to the high internal precision of the atmospheric pa¬ 
rameters derived in our solar twin stars, we were able to employ 
standard isochrone techniques to measure reliable relative ages 
for these objects. We adopted the ages previously derived by our 
group (Monroe et al. 2013, 2015; Melendez et al. 2014a) using 
the algorithm described in Ramirez et al. (2013, 2014), which 
computes the age probability distribution function by compar¬ 
ing the location of the star on the r e ff,logg, [Fe/H] parameter 
space with the values predicted by theory. For the two youngest 
stars, we also used other age indicators, as described in Mon¬ 
roe et al. (2015). In short, we adopted a rotational age (0.5+0.2 
Gyr) for HD 20630, which excellently agrees with iRibas et all 
d2010l) . who obtained 0.6+0.2 Gyr using different indicators. For 
HD 202826 a rotation period is not available to estimate its age, 
hence we used its chromospheric activity. X-ray luminosity, and 
isochrones. From the isochrones we derive an average mass for 
the sample stars of 1.01 ± 0.03 M Q , that is, solar within the er¬ 
rors, reinforcing thus the use of solar twins as proxies of the Sun 
at different ages. The estimated ages are given in Table 1.0 
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Fig. 2. Be (top) and Li (bottom) abundances vs. age. For the mod¬ 
els of Be depletion, we adopt an i nitial m eteoritic A(Be) = 1.45 dex, 
which is based on the value bv iLodders ( 2003 ) (1.41 dex), plus the 
0.04 dex gravitational settling effect over 4.5 Gyr. The models of Li 
depletion were normali zed to the solar Li abunda nce. The green solid 
lines are the models by Pinsonneault et alj dl989l). the red dotted lines 
represent the models bv Ido Nascjmentneta a (2009), the blue dashed 
lines are the modified models of ldo Nasciinentoet all d2009l) (see text), 
and the purple long-d ashed lines are pre dictions of Be depleti on by 
IXiong & De ng (2007) and Li depletion bv IXiong & Penal (2009). The 
teal dotted-dashed line (top panel) is the same Be depletion model of 
Ido Nascimento et ahl d2009l) shown by the red dotted lines, but with a 
higher initial A(Be)= 1.50 dex (rather than A(Be) = 1.45 dex). The ad¬ 
ditional 0.05 dex is to compensate for the refractory depletion of Be in 
the Sun (see text). 


4. Discussion 

Figure [2] (top panel) shows the measured Be abundances in our 
solar twins as a function of stellar age. The scatter in the Be con¬ 
tent is very small: 0.04 dex. This is significantly less than in pre¬ 
vious studies; for example. [Takeda et alJd2011l) analyzed Be in a 
large sample of solar analogs and found a dispersion of 0.2 dex. 
A linear fit to our data using the error bars in both age and Be 
abundances gives a slope of -8.09 • 10~ 3 ±4.17 ■ 10 -3 dex/Gyr, 
which is a shallow trend at best. For comparison, we show in 
the bottom panel the corresponding NLTE Li abundances for the 
same sample stars (Monroe et al. 2015). A fit to the Li data gives 
a slope of -0.23+0.01 dex/Gyr. Clearly, there is a steep Li deple¬ 
tion, but the mixing processes responsible for the destruction of 
lithium must be relatively shallow and cannot transport material 
to deeper regions where significant Be destruction can occur. 

In standard evolution models of the Sun, the depletion of 
the light elements Be and Li is expected to occur only below 

2 Note that our ages are differential and that our errors are only inter¬ 
nal. We caution that the ages may turn out to be slightly older or younger 

depending on the choice of isochrones, making the age range covered by 
the sample slightly wider or narrower, but the relative ages are reason¬ 
ably well constrained, as shown for example in Fig. 5 of lMelendez et all 
d2014al) , where the relative ages between the Sun and the solar twin 18 
Sco are consistent for two different sets of isochrones. 


the convective zone and thus the surface abundances should re¬ 
main unchanged, obviously in contrast to the observational evi¬ 
dence in the case of Li. Only for less massive stars do standard 
stellar models predict significant Be depletion as a result of the 
deeper convective zone (Sant os et al.l2004l). Owing to the mono - 
tonically increasing Li depletion with age dMonroe et al.ll2013l) . 
most Li destruction occurs on secular timescales during the main 
sequence. Since standard models cannot predict this behavior, 
there must be extra mixing below the convection zone that brings 
material down to sufficiently large depths and temperatures, and 
then back into the convection zone. In Fig. [2] four different pre- 
dictions for Be and Li d epletion are compared. The model by 
[Pinsonneault et al] ( 1989 ) considers rotationally induced mixing 
and depletes Li reasonably wel l, but depletes far too much Be 
(green solid line). The model of ldo Nascimento et ali ( 2009 ) in¬ 
cludes extra mixing due to diffusion (including gravitational set¬ 
tling) and rotation. It reproduces Li in solar twins well, but Be 
depletion is somewhat larger than observed (red dotted lines). 
We have modified this model to calibrate the amount of merid¬ 
ional circulation, achieving a steep and a shallow Li and Be 
depletion, respectively (blue dashed lines). The new model has 
meridional circulation with a lower efficiency, hence destroying 
less Be and Li, and we increased the turbulent diffusion coeffi¬ 
cient of the tachocline below the convective zone to destroy more 
Li (without affecting Be). Finally, the models of lXiong & Pend 
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( 2007 . 120091) incorporate convective overshoot as well as grav¬ 
itational settling. They reproduce the Li depletion with age but 
do not deplete Be (purple long-dashed line), unlike the observa¬ 
tions, which seems to suggest a shallow depletion of Be. 

The initial (zero-age) Be abundance for the model predic¬ 
tions assumes them to be equal to the meteoritic values. For the 
meteoritic abundance, iLoddersI ( 200 3) recommended A(Be)= 
1.41 + 0.08 dex, which was su bseque ntly revised to A(Be)= 
1.32 + 0.03 dex bv lLodders et al.l (120091) . These are indirect mea¬ 
surements, however, because Be is difficult to measure in car¬ 
bonaceous chondrites of type Cl, which are the least modified 
meteorites and thus the preferred choice when inferring the pri¬ 
mordial solar system abundances. Instead, the meteoritic value 
was estimated from the relative abundances of refractory ele¬ 
ments in CM and CV chondrites in which Be has been measured 
in a couple of cases. The relatively large uncertainty for Be for 
being meteorites reflects this indirect procedure. 

Since the absolute abundance scales for meteorites are set 
by enforci ng that the photospheric and meteori tic Si abundances 
are equal (Asplund 200(1 lAsplund et aPl2009h and all elements 
heavier than hydrogen have experienced gravitational settling in 
the Sun over the past 4.5 Gyr (see discussion in lAsnlund et alj 
l2009h , a more appropriate initial Be value for these model pre¬ 
dictions would thus be 0.04 de x hi gher than the abund ances rec¬ 
ommended bv lLoddersl (120031) and Lodders et alJ (l2009h . Our so¬ 
lar twins data would seem to suggest that the higher meteoritic 
Be abundance (1.41 dex) is more appropriate. We thus added 
0.04 dex to this value, adopting A(Be) = 1.45 dex as the initial 
Be abundance in Fig. [2] The Sun may be slightly less abundant 
in Be by ~ 0.05 dex for its age compared with other similar solar 
twins (Fig. 0, but a larger number of solar twins would be re¬ 
quired to confirm this impression. Perhaps the somewhat lower 
solar Be abu ndance could arise sim ply because the Sun is poor in 
refractories (IMelendez et al. l2009h: Be has a condensation tem¬ 
perature of 7' conl i = 1452 K ( Loddersll2003l) . which means that it 
is a refractory element. The Sun is probably deficient in refrac¬ 
tories as a result of the formation of rocky planets in the solar 
system (iMelendez et al.|[2009t) . For its 7 a)ncl , we estimate that Be 
should be depleted in ~0.05 dex in the Sun. Interestingly, if we 
consider an initial A(Be) =1.5 dex owing to the refractory deple¬ 
tion of Be, the modified model by Ido Nascimento et al] ( i2009i) 
reproduces the Be abundances of most solar twins well (dotted- 
dashed line in Fig.0. 

Recently, lAdibekvan et alJ d2014l) suggested that the deple¬ 
tion of refractory elements in the Sun relative to solar twins could 
be an age effect. If this interpretation is correct, and because Be 
is a refractory element, we should have found that beryllium in 
solar twins older than the Sun (4.6 Gyr) is depleted relative to 
the Sun, because lAdibekvan et al.l (12014 ) analyzed their solar 
analogs relative to the Sun. However, this is not what we ob¬ 
serve in Fig. 0 Thus, our results seem to be in conflict with the 
interpretation by lAdibekvan et akl ( 2014 ). 


5. Conclusions 

We presented the first detailed study of beryllium abundances 
in solar twins covering a broad range of ages (0.5 - 8.2 Gyr). 
Our analysis revealed that the Be abundance is relatively con¬ 
stant with age, with a scatter of only 0.04 dex and a weak, if any, 
trend with age. This is in contrast to the large observed deple¬ 
tion of Li with age (Baumann et al.l 120101 : iMonroe et al .11201 31 : 
IMelendez et alJl2014bh . showing that the transport mechanisms 
are deep enough to reach the region where Li is burned, but not 
deep enough to reach the higher temperatures needed to burn 


Table 1. Be abundances and errors for the eight solar twins and the Sun, 
together with the inferred macroturbulence, v sin i, and ages of the stars. 


Star 

^macroA 7 $in i 

A(Be) 

param a /ob s b /total c 

Age 


(km s -1 ) 

(dex) 

(dex) 

(Gyr) 

HD20630 

3.5/4.2 

1.44 

0.03/0.03/0.04 

0 5 +u ^ 

0 9 + ° b 
-Q-5 

6 0 + °- 5 

HD202628 

3.7/2.4 

1.51 

0.01/0.04/0.04 

HIP30502 

3.5/1.6 

1.43 

0.01/0.03/0.03 

HIP73815 

3.6/1.7 

1.49 

0.01/0.03/0.03 

HIP77883 

3.4/1.8 

1.42 

0.01/0.02/0.02 

2 9 +I * 
«%? 

HIP89650 

3.8/1.7 

1.47 

0.01/0.02/0.02 

18Sco 

3.7/2.0 

1.52 

0.01/0.03/0.03 

HIP102152 

3.5/1.8 

1.41 

0.01/0.04/0.04 

Sun 

3.6/1.9 

1.38 

0.00/0.01/0.01 

4.6 


a errors due to stellar parameters 
b observational errors 

c quadric sum of the observational and systematic errors 


Be. Our Li and Be results provide stringent constraints on stellar 
models and nonstandard mixing processes beyond treating con¬ 
vection through the mixing length theory. 
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